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Abstract — The present study deals with the combined effects 
of Soret and chemical reaction on the radiative MHD free 
convection flow of an electrically conducting incompressible 
viscous fluid over an inclined plate embedded in a porous 
medium is studied. The impulsively started plate with variable 
temperature and mass diffusion is considered. The 
dimensionless governing partial differential equations have been 
solved numerically and expressions for velocity, temperature 
and concentration profile are obtained. They satisfy all imposed 
initial and boundary conditions and can be reduced, as special 
cases, to some known solutions from the literature. Expressions 
for skin friction, Nusselt number and Sherwood number are also 
obtained. The effects of pertinent parameters on velocity, 
temperature and concentration profiles are graphically 
displayed where as the variations in skin friction, Nusselt 
number and Sherwood number are shown through tables. 


Index Terms — Soret number (Sr), Inclined plate, MHD free 
convective, Viscous flow, Chemical reaction 


I. Introduction 

The conjugate effects of heat and mass transfer occurs as a 
result of combined buoyancy effects of thermal and mass 
diffusion which essentially plays an important role in 
Geophysics, Aeronautics and Chemical Engineering. Guo Z, 
Sung Hj [1] studied the conjugate effect of heat and mass 
transfer in metal hydride beds in the hydriding process which 
is useful in food drying, food processing and polymer 
production. The conjugate phenomenon of heat and mass 
transfer either analytically or numerically gained a 
considerable amount of interest in recent years and extensive 
research work is done on account of this [2-5]. Bhuvaneswari 
M et al. [6] studied an exact analysis of radiation convective 
flow heat and mass transfer over an inclined plate in a porous 
medium. Kandasamy R et al. [7] analysed the Lie group 
analysis for the effect of temperature-dependent fluid 
viscosity with thermophoresis and chemical reaction on MHD 
free convective heat and mass transfer over a porous 
stretching surface in the presence of heat source/sink. 

Ahmed N and Sarmah H.K [8] investigated the thermal 
radiation effect on a transient MHD flow with mass transfer 
past an impulsively fixed infinite vertical plate. Rajput and 
Kumar [9] thoroughly studied the MHD flow past an 
impulsively started vertical plate with variable temperature 
and mass diffusion. They used the Laplace transform method 
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to find the exact solutions for velocity, temperature and 
concentration. In a subsequent year Rajput and Kumar [10] 
extended by taking the thermal radiation effect. Ahmad [11] 
studied MHD transient free convection and mass transfer flow 
of a viscous, incompressible and electrically conducting fluid 
in the presence of thermal diffusion and thermal radiation. He 
obtained exact solutions for velocity, temperature, and 
concentration using the Laplace transforms method. Osman et 
al. [12] studied analytically the thermal radiation and 
chemical reaction effects on unsteadyMHD free convection 
flow in a porous medium with heat source/ sink. 

By taking the porous medium effect, Sami et al. [13] 
provided an exact analysis to the study of the 
magnetohydrodynamic free convection flow of an 
incompressible viscous fluid past an infinite vertical 
oscillating plate with uniform heat flux. In other investigation, 
Sami et al. [14] studied the MHD free convection flow in a 
porous medium with thermal diffusion and ramped wall 
temperature. They obtained exact dimensionless solutions of 
momentum and energy equations, under Boussinesq 
approximation using the Laplace transforms. Magyari E, 
Pantokratoras A [15] gave a note on the effect of thermal 
radiation in the linearized Rosseland approximation on the 
heat transfer characteristics of various boundary layer flows. 

Ziyauddin and Kumar [16] studied the radiation effects on 
unsteady MHD natural convection flow in a porous medium 
with conjugate heat and mass transfer past a moving inclined 
plate in the presence of chemical reaction, variable 
temperature and mass diffusion. They used an explicit finite 
difference method to solve the coupled linear partial 
differential equations numerically, and the results are 
graphically displayed. Unfortunately, in this work the plate is 
not porous as the authors mentioned in the paper. On the other 
hand the numerical solutions of the free convection problems 
are more convenient and easy to handle as compare to exact 
solutions. Muthucumaraswamy and Janakiraman[17] 
obtained exact solutions for the MHD flow of viscous 
optically thin fluid past a vertical flat plate in a non-porous 
medium. They considered uniform heat and variable mass 
transfer. 

Larhad Ali [18] presented an exact analysis of combined 
effects of radiation and chemical reaction on the MHD free 
convection flow of an electrically conducting incompressible 
viscous fluid over an inclined plate embedded in a porous 
medium. The impulsively started plate with variable 
temperature and mass diffusion is considered. The 
dimensionless momentum equation coupled with the energy 
and mass diffusion equations are analytically solved using the 
Laplace transform method. 

The aim of this study is to present an analysis of combined 
effects of Soret and viscosity the radiative MHD free 
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convection flow of an electrically conducting incompressible 
viscous fluid over an inclined plate embedded in a porous 
medium. The impulsively started plate with variable 
temperature and mass diffusion is considered. The 
dimensionless momentum equation coupled with the energy 
and mass diffusion equations are numerically solved and 
expressions for velocity, temperature and concentration fields 
are obtained. They satisfy all imposed initial and boundary 
conditions and can be reduced, as special cases, to some 
known solutions from the literature. Expressions for skin 
friction, Nusselt number and Sherwood number are also 
obtained. Finally, the effects of pertinent parameters on 
velocity, temperature and concentration profiles are 
graphically displayed whereas the variations in skin friction, 
Nusselt number and Sherwood number are presented clearly 
through tables. 

II. Formulation of the Problem 

A radiative MHD free convection flow of an electrically 
conducting incompressible viscous fluid over an inclined 
plate embedded in a porous medium is considered under the 
following assumptions 

1 . The plate is started impulsively with variable 

temperature and the mass diffusion is considered 

2. The x - axis is taken along the plate with the angle of 

inclination (a) to the vertical and the y - axis is taken 
normal to the plate. 

3. The viscous fluid is taken to be electrically conducting 

and fills the porous half space y > 0. 

4. A uniform magnetic field of strength B 0 is applied in 

the y -direction transversely to the plate. The applied 
magnetic field is assumed to be strong enough so that 
the induced magnetic field due to the fluid motion is 
weak and can be neglected. 

5. According to Cramer and Pai[5], this assumption is 

physically justified for partially ionized fluids and 
metallic liquids because of their small magnetic 
Reynolds number. Since there is no applied or 
polarization voltage imposed on the flow field, the 
electric field due to polarization of charges is zero. 


6. Initially, both the fluid and the plate are at rest with 

* 

constant temperature T x and constant 

concentration C rf , . At time t = 0 , the plate is 

given a sudden jerk, and the motion is induced in the 
direction of flow against the gravity with uniform 
velocity u c .The temperature and concentration of 
the plate is raised linearly with respect to time. 

7. The fluid is thick gray absorbing-emitting radiation 

but non-scattering medium. 

8. Since the plate is infinite in the (x , z ) plane, all 

physical variables are functions of y and t* only. 

The physical model and the coordinate system is shown in the 
following figure 
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In view of the above assumptions, as well as of the usual Boussineq’s approximation, the governing equations reduce 
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Following Magyari and Pantokratoras [15] the Rosseland’s approximation for radiative flux ( tj r ) is adopted 
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The temperature differences within the flow are assumed to be so small and T can be expressed as a linear function of the 
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Into the governing equations and initial and boundary conditions to get the governing equations with the corresponding initial 
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used to solve the governing equations subject to the respective 
III. Finite Difference Scheme initial and boundary conditions. Using the finite difference 

Since exact or approximate solutions are not possible for scheme ’ the governing equations with the corresponding 
this set of equations, an implicit finite difference method is boundary conditions become 
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The suffixes, i corresponds to y and j corresponds to t and At = / ( 7 + 1 ) — t (i) and Ay = y(j + l) ~ y(j). The 
computations were carried out for different values of the various physical parameters. 


The Skin - friction coefficient or the shear stress at the wall is given by 

du(y, t ) 

T = - 

6y 

L y = 0 

The rate of heat transfer in terms of the Nusselt number is 
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The rate of mass transfer in terms of Sherwood number is 

dt/)(y, t ) 


Sh = 


dy 


(19) 


(20) 


( 21 ) 


IV. Discussion of Results 


y = 0 

system parameters. The results are shown in graphs and 
tables. 

The velocity profiles are discussed through graphs 1 to 6 for 

various parameters such as Hartmann number (M), thermal 
The numerical solutions for the velocity, temperature, and Grashoff number ( G r), the Solutal 
mass diffusion are computed for various physical parameters 
such as Hartmann number (M), Eckert number (Ec), Prandtl 
number (Pr), time (7), and Schmidt number (Sc), etc. The skin 
- friction, the rate of heat transfer in terms of the Nusselt 
number (Nu), and the rate of mass transfer in terms of the 
Sherwood number (Sh) are also derived in terms of the given 
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Fig.l 

Variation of velocity profiles for different values of M 



Fig. 5 

Variation of velocity profiles for different values of Sc, y and a 



Fig. 2 Velocity profiles for different values of Gr 



Fig. 3 

Variation of velocity profiles for different values of Gr and Go 



Fig. 4 

Variation of velocity profiles for different values of Pr and K 



Fig. 6 

Variation of velocity profiles with inclination a 

The effect of Magnetic parameter M on the velocity profile is 
presented in figure 1 . It is seen that the velocity decreases as 
the magnetic parameter M is increases. Physically, it is 
justified because the application of transverse magnetic field 
always results in a resistive type force called Lorentz force 
which is similar to drag force and tends to resist the fluid 
motion, finally reducing its velocity. 

The variation of velocity distribution to the thermal 
Grashoff number is discussed in the figure 2. The thermal 
Grashoff number is the relative effect of the thermal buoyancy 
force to the viscous hydrodynamic force in the boundary 
layer. The positive values of Grashoff number indicates the 
cooling of the plate. The rise in the velocity is observed due to 
the enhancement of the thermal buoyancy force. As thermal 
buoyancy increases, the velocity increases rapidly near the 
plate and gradually decreases to free stream velocity. With the 
decrease of the thermal Grashof number i.e Gr<0, from figure, 
an opposite behavior is observed. 

The combined effect of Grashof number (Gr) and the 
Solutal Grashof number (Gc) is shown in the figure 3. 

Similar to the thermal Grashoff number, the 
Solutal Grashoff number effect is also to increase in the 
velocity. The rise in the velocity distribution with the increase 
of Solutal Grashof number observed in the figure and also the 
gradual decrease in the boundary layer. 

In the figure 4, the combined effect of Prandtl number (Pr) 
and the permeability parameter (K) is given. Physically, it is 
true due to the fact that an increase in Prandtl number 
increases the viscosity of the fluid, becomes thick and 
consequently leads to a decrease in the velocity. But increase 
in permeability increases the size of the pores which 
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influences the flow to be fast so that the velocity increases 
with the increase of Pr and K. 

The variation of velocity distribution with the physical 
parameters Schmidt number (Sc) and the chemical reaction 
parameter (y) is observed from the figure 5. With the increase 
of Sc and y, the velocity rises but decreases gradually in the 
boundary layer and reaches a steady state. This effect can be 
seen from the figure. It is noticed from the figure 6 that the 
velocity increases sharply as the inclination decreases. 



Fig. 7 

Variation of temperature profiles for different values of 
Pr and K 



Y 


Fig. 8 

Variation of temperature profiles for different values of 
Pr and N when t=2 



Y 


Fig. 9 

Variation of temperature profiles for different values of t 
when Pr=0.71 and N=1 



Y - 

Fig. 10 

Variation of temperature profiles for different values of t 
when Sc=2 and y=l 



Y 


Fig. 1 1 

Variation of temperature profiles for different values of Ec 

Figure 7 depicts the influence of Prandtl number (Pr) and 
permeability parameter (K) on the temperature. Four different 
values of Prandtl number namely Pr~0:71, 1, 7, and 11 
corresponding to air, electrolyte, water and lubricants 
respectively are chosen. It is mainly observed that 
temperature decreases with increasing Pr. It is interesting to 
note that a comparison of curves shows that the temperature 
falls more rapidly for water in comparison to air. A similar 
behaviour is also expected due to the fact that increasing Pr 
means the increasing viscous nature of a fluid. For smaller 
values of Prandtl number fluids possess high thermal 
conductivity and heat diffuses away from the surface faster 
than at higher values of Prandtl number. Thus the boundary 
layer becomes thicker and consequently the temperature 
decreases when Pr is increased. Furthermore, the temperature 
profiles for increasing values of permeability parameter K 
indicate a decreasing behaviour as shown in the figure. 

The temperature profiles for increasing values of radiation 
parameter (N) and Prandtl number (Pr) is shown in the figure 
8. It is expected to observe an increase in the temperature but 
with the simultaneous increase in the Prandtl number 
decreases the temperature. This can be seen from the figure. 

Figures 9 and 10 is plotted to show the effect of the 
dimensionless time t on the temperature when 
Pr = 0.7 1 and /V = I and Sc = 2 and ;/ = I Four 

different values of time t =0.4, 0.6, 0.8, and 1.0 are 
chosen. Obviously, the temperature increases with increasing 
time. 

The rise of the temperature distribution with the Eckert 
number is shown in the figure 1 1 . The temperature increases 
with the respective increase in the Eckert number. The 
temperature increases slightly near the plate but decreases 
sharply away from the plate i.e. with the increase in the 
boundary 
layer. 



Fig. 12 

Variation of Concentration profiles for different values of Sr and y 
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Fig. 13 

Variation of Concentration profiles for different values of Pr and K 



Fig. 14 

Variation of Concentration profiles for different values of 
Sc, y and a 



Y 


Fig. 15.1 

Variation of Concentration profiles for different values of t when 
Pr=0.71 and N=1 



Fig. 15.2 

Variation of Concentration profiles for different values of t 
when Sc=2 and y=l 



Fig. 17 

Variation of Concentration profiles for different values of Sr 



Fig. 18 

Variation of Concentration profiles for different values of Ec 

Figure 12 shows the influence of the Soret number (Sr) and 
chemical reaction parameter (y) on the concentration profiles. 
The rise in the concentration is observed here for the increase 
of Soret number (Sr) and chemical reaction parameter (y). 

In the figure 13, the variation of concentration, with the 
simultaneous increase of the Prandtl number (Pr) and 
permeability parameter (K), is shown. The effect of 
permeability parameter is to decrease the concentration since 
with the increase of permeability, the fluid becomes thinner 
but with the simultaneous increase of Prandtl number, the 
concentration increases. Prandtl number increases the 
viscosity which in turn increases the concentration. 

Graphical results of concentration profiles for different 
values of Schmidt number Sc and chemical reaction 
parameter (y) are shown in Figure 14. The figure shows that an 
increase in Sc, decreases the concentration. Further more, it is 
interesting to note that the concentration profiles fall slowly 
and steadily for Hydrogen (Sc~0:2) and Helium (Sc~0:3) but 
falls very rapidly for water vapour (Sc~0:6). Physically this is 
true because of the fact that the water vapour can be used for 
maintaining normal concentration field where as Hydrogen 
can be used for maintaining effective concentration field. The 
physical effect of chemical reaction parameter (y) clearly 
demonstrates that concentration profiles decrease rapidly 
when y is increased. 

Figures 15.1 and 15.2are plotted to show the comparative 
study of the dimensionless time t on the concentration when 
Pr = 0.7 1 and N — l and Sc = 2 and y = 1 Four different 
values of time t = 0.4, 0.6, 0.8, and 1.0 are chosen. We 
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found that the concentration increases when the time (t) is 
increased. 

The variation of concentration distribution with the 
physical parameters Schmidt number (Sc) and the chemical 
reaction parameter (y) is observed from the figure 16. With 
the increase of Sc and y, the velocity decreases gradually in 
the boundary layer and reaches a steady state. This effect can 
be seen from the figure. 

The effects of Soret and Eckert numbers on the 
concentration distribution are plotted in the figures 17 and 18. 
As the Soret number increases the concentration increases in 
the middle of the boundary layer but eventually decreases. 
This can be observed from the figure 17. The gradual thinning 
of the fluid in the boundary layer with the increase of Eckert 
number (Ec) is showed in the figure 18. 


Table-I Variation in Skin-friction coefficient(x) 


M 

r 

Gr 

Gc 

Pr 

K 

Sc 

a 

x(Skin-f 

riction) 

5 

i 

l 

1 

0.71 

0.1 

0.22 

-n/b 

- 2.273828 

10 

i 

l 

1 

0.71 

0.1 

0.22 

-n/b 

- 3.126577 

5 

3 

1 

1 

0.71 

0.1 

0.22 

-nib 

- 2.491116 

5 

1 

3 

1 

0.71 

0.1 

0.22 

-n/b 

- 1.853089 

5 

1 

1 

4 

0.71 

0.1 

0.22 

-nib 

- 1.958154 

5 

1 

1 

1 

1 

0.1 

0.22 

-n/b 

- 2.384553 

5 

1 

1 

1 

0.71 

0.5 

0.22 

-nib 

- 2.672829 

5 

1 

1 

1 

0.71 

0.1 

0.60 

-n/b 

- 2.491116 

5 

1 

1 

1 

0.71 

0.1 

0.22 

n/b 

- 1.60452 


The skin-friction coefficient at the wall (x) is given in the 
tables I. It is clearly noticed that the increase in magnetic field 
(M), permeability parameter (K), Solutal Grashof number 
(Gc) and Prandtl number (Pr) causes the decrease in 
Skin-friction coefficient. The shear stress increases with the 
increasing values of the chemical reaction parameter (y), 
thermal Grashof number (Gr), Schmidt number (Sc) and the 
inclination of the plate(a) . But for constant values 
ofEc.Nand.t. 


Table II Variation of rate of heat transfer (Nu) 


M 

Y 

Gr 

Gc 

Pr 

K 

Sc 

N 

t 

Ec 

a 

Nu 

5 

1 

1 

1 

0.71 

0.1 

0.22 

1 

0.4 

1 

- n /6 

-.8215216 

10 

1 

1 

1 

0.71 

0.1 

0.22 

1 

0.4 

1 

- n /6 

-.7780302 

5 

3 

1 

1 

0.71 

0.1 

0.22 

1 

0.4 

1 

- n /6 

-.480635 

5 

1 

3 

1 

0.71 

0.1 

0.22 

1 

0.4 

1 

- n /6 

-.9846448 

5 

1 

1 

4 

0.71 

0.1 

0.22 

1 

0.4 

1 

- n /6 

-.4840673 

5 

1 

1 

1 

1 

0.1 

0.22 

1 

0.4 

1 

- n /6 

-.4738897 

5 

1 

1 

1 

0.71 

0.5 

0.22 

1 

0.4 

1 

- n /6 

-.471469 

5 

1 

1 

1 

0.71 

0.1 

0.60 

1 

0.4 

1 

- n /6 

-.9806249 

5 

1 

1 

1 

0.71 

0.1 

0.22 

2 

0.4 

1 

- n /6 

- 2.404324 

5 

1 

1 

1 

0.71 

0.1 

0.22 

1 

0.6 

1 

- n /6 

-.5711901 

5 

1 

1 

1 

0.71 

0.1 

0.22 

1 

0.4 

3 

— n /6 

-.418885 

5 

1 

1 

1 

0.71 

0.1 

0.22 

1 

0.4 

1 

n /6 

-.8062888 


The table III gives the variation of rate of heat transfer in 
terms of Nusselt number. With the increase of magnetic field 
(M), Solutal Grashof number(Gc), Chemical reaction 
parameter (y), Prandtl number (Pr), Permeability 
parameter(K),time t and Eckert number (Ec) causes the 


increase in the Nusselt number. The Nusselt number 
decreases with increase in thermal Grashof number (Gr), 
Schmidt number (Sc), N and inclination of the plate (a). 


Table III Variation of rate of mass transfer (Sh) 


M 

y 

Gr 

Gc 

Pr 

K 

Sc 

N 

t 

Ec 

i 

Sh(Sher.vood 

number) 

5 

l 

1 

1 

0.71 

0.1 

0.22 

1 

0.4 

1 

-ir/6 

-1.981967 

10 

l 

1 

1 

0.71 

0.1 

0.22 

1 

0.4 

1 

-jt/6 

-2.025409 

5 

3 

1 

1 

0.71 

0.1 

0.22 


0.4 

1 

-n/6 

-.6178092 

5 

1 

3 

1 

0.71 

0.1 

022 


0.4 

1 

-n/6 

-1.073927 

5 

1 

i 

4 

0.71 

0.1 

022 

1 

0.4 

1 

-n/6 

-.6082868 

5 

1 

i 

1 

1 

0.1 

0.22 

1 

0.4 

1 

-n/6 

-.6122832 

5 

I 

i 

1 

0.71 

0.5 

0.22 

: 

0.4 

: 

-n/6 

,5392113 

5 

1 

i 

1 

0.71 

0.1 

0.60 

l 

0.4 

l 

■n/6 

-1.381791 

5 

1 

i 

1 

0.71 

0.1 

0.22 

L 

0.4 

l 

-n/6 

-1.973262 

5 

I 

i 

1 

0.71 

0.1 

0.22 

1 

0.6 

l 

■n/6 

-1.114278 

5 

I 

i 

1 

0.71 

0.1 

0.22 

1 

0.4 

3 

-n/6 

-2.332864 

5 

1 

i 

1 

0.71 

0.1 

0.22 

1 

0.4 

i 

n/6 

-1.997039 


From the table III, it is observed that with the increase of the 
magnetic field (M), Eckert number (Ec), thermal Grashof 
number (Gr), Prandtl number (Pr)and Schmidt number (Sc) 
leads to decrease in Sherwood number, the inclination of the 
plate, the rate of mass transfer increases with the increase of 
chemical reaction parameter (y), Solutal Grashof number 
(Gc), Permeability parameter (K), time t and inclination of the 
plate (a). 

Particular cases 

The following published results are reduced as special cases 
from the present solutions. 

1. The results obtained in the absence of viscosity, 

thermal and mass transfer effects are quite 
identical to the known result obtained 
bySamiulhaqet.al [13] for CO —> 0 the 
impulsive motion of the plate. Furthermore, the 
solution corresponding to hydrodynamic fluid 
passing through a non - porous medium is 
obtained by neglecting the magnetic and porous 
effects. 

2. In the absence of variable mass diffusion and when 

Sr = 0 the results are similar with the results of 
Muthucumaraswamy R, Janakiraman B[17]. 

3. In the absence of radiation, viscous and mass 

diffusion and, the solutions through the 
non-porous medium, obtained were in agreement 
with the results of Rajput and Kumar (9). 

V. Conclusions 

An analysis of combined effects of Soret and Chemical 
reaction on the radiative MHD free convective flow of an 
electrically conducting incompressible viscous fluid over an 
inclined plate embedded in a porous medium is considered. 
The dimensionless momentum equation coupled with the 
energy and mass diffusion equations are numerically solved. 
The effects of pertinent parameters on velocity, temperature 
and concentration are graphically displayed where as the 
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variations in skin friction, Nusselt number and Sherwood 
number are presented through tables, we conclude that 

1 . The velocity profiles increases with increase in Gr, Gr 

&Gc, Pr &K, Sc &y and decrease with increase in M, 
Gr<0 and a. 

2. Temperature profiles increases with increase in t & Ec 

and decrease with increase in Pr& K, and Pr& N. 

3. Increase in Sr&y,Sc&y,Ec leads to decrease in 

concentration profiles and increase in Sr, Sc &t 
causes increase in concentration profiles. 

4. With an increase in M, Gc, Pr and K tends to decrease 

in Skin-friction while increase in y, Gr, aSc and 
leads to increase in Skin-friction. 

5. Increase in M, Gc, Pr, K,t, Ec and y causes decrease in 

Nusselt number and Increase in Gr, Sc, N and a 
leads to decrease in Nusselt number. 

6. The Sherwood number increases with increase in Gc, 

K, t, y and a while decease with increase in M, Gr, 
Pr, Sc, N and Ec. 

7. The results obtained in the absence of viscosity, 

thermal and mass transfer effects are quite identical 
to the known result obtained by Samiulhaq et al. [13] 
lor 0) — > 0 the impulsive motion of the plate. 
Furthermore, the solution corresponding to 
hydrodynamic fluid passing through a non - porous 
medium is obtained by neglecting the magnetic and 
porous effects. 

8. In the absence of variable mass diffusion and 

when Sr = 0 the results are coincides with the 
results of Muthucumaraswamy R, lanakiraman 
B [ 17] . 

9. In the absence of radiation, viscous and mass diffusion 

and, the solutions through the non - porous medium, 
obtained were in agreement with the results of 
Rajput and Kumar [9], 
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